Evolving network models under a dynamic growth rule 
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Evolving network models under a dynamic growth rule which comprises the addition and deletion 
of nodes are investigated. By adding a node with a probability P a or deleting a node with the 
probability Pd = 1 — P a at each time step, where P a and Pd are determined by the Logistic population 
equation, topological properties of networks are studied. All the fat-tailed degree distributions 
observed in real systems are obtained, giving the evidence that the mechanism of addition and 
deletion can lead to the diversity of degree distribution of real systems. Moreover, it is found 
that the networks exhibit nonstationary degree distributions, changing from the power-law to the 
exponential one or from the exponential to the Gaussian one. These results can be expected to shed 
some light on the formation and evolution of real complex real- world networks. 
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I. INTRODUCTION 

Considerable interest is focused on complex networks 
currently due to their potential to describe many systems 
in nature and society [l|, [2| . Theoretical models are de- 
veloped to reproduce topological properties of these sys- 
tems The recently observed scale-free (SF) prop- 
erty Q has made the scientists to realize that static net- 
work models 0, Q do not provide appropriate descrip- 
tions to real systems which essentially keep growing with 
time [J, Q . Barabasi and Albert [5j proposed a simple 
evolving network model (BA model) to explain such SF 
property. In this model, the growing nature of real sys- 
tems is captured by a BA-type growth rule. According 
to this rule, one node is added into the network at each 
time step, intending to mimic the growing process of real 
systems. Another ingredient in this model is the mecha- 
nism of preferential attachment (PA) , which assumes that 
newly added nodes are attached preferentially to nodes 
with higher degrees @. Based on the BA-type growth 
rule, many evolving network models are introduced 0, 0] ■ 
The modeling framework of network evolution is formed. 

Among the many quantities proposed to characterize 
topological properties of networks 0, 0], the degree dis- 
tribution p(k), which gives the probability that a node 
in the network possesses k edges, is of particular impor- 
tance. The study of networks has undergone a transi- 
tion from the investigation of static network models with 
Poisson degree distribution 0, Q to the explanation of 
real systems with various fat-tailed degree distributions 
(FTDDs) [Ji]. It was found that many real- world net- 
works are characterized by the power-law degree distri- 
bution (PLDD) [ah as well as the exponential degree dis- 
tribution (EDD) 7]. Some more exhaustive experiment 
reveals that, in addition to the PLDD and the EDD, 
the truncated power-law degree distribution (TPLDD) 
0, Q and the truncated exponential degree distribution 
(TEDD) [T(J are also observed. Furthermore, in Ref. [8J, 
the Gaussian degree distribution (GDD) was reported as 
well. 



One theoretical challenge has been to explain the origin 
of these observed FTDDs. In the BA model, the fact is 
revealed that growing networks with PA yield the PLDD 
In addition, growing networks without PA were also 
studied in the random evolving network model (REN 
model) [ll|, [HI , in which networks grow by the growth 
rule of BA-type, while the newly added nodes connect to 
randomly chosen existing ones. In this model, the EDD 
was obtained. In order to explain those observed degree 
distributions which are neither strict power-law nor strict 
exponential [81-4101] . many other mechanisms were added 
into this two models, such as the mechanism of adding 
and rewiring edges between existing nodes [IBJIrjj . the 
mechanism of ageing and cost Q, as well as the mecha- 
nism of information filtering Q . These more specialized 
models indeed created the TPLDD or the EDD in some 
parameter regimes. However, relatively lesser attention 
is paid to systematical studies on the origin of the diver- 
sit y of the observed FTDDs, with some recent exceptions 

The BA-type growth rule gives a somewhat simplified 
description to the evolution of real system. As a matter 
of fact, in real growing networks, there are constant ad- 
dition of new elements, but accompanied b y p ermanent 
removal of old elements (deletion of nodes) [l9l - l2l| . The 
scaling behavior of growing networks has already been 
shown to be strongly affected when the deletion of node 
is taken into account [22U25| . In a recent work [13], we 
proposed a new type of network growth rule which com- 
prises the addition and deletion (AD) of nodes. Based 
on such AD growth rule, by adding a node with a proba- 
bility P a or deleting a node with probability Pd = 1 — P a 
at each time step, topological properties of growing net- 
works with and without PA are studied, respectively. In 
ref. [13], we considered a simple case: P a and Pd were 
treated as constants. [P a is an adjustable parameter in 
the model.) All the observed FTDDs of real systems were 
obtained in the model, indicating that the mechanism of 
AD can lead to the diversity of FTDD in real systems. 

The AD growth rule introduced in [ljj gives a more de- 
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tailed description to the evolution of real systems. How- 
ever, this evolution can be even more complex. For ex- 
ample, as real systems grow, limited resources will cause 
interelement competition which, in turn, has a tendency 
to retard the growth of these systems. Such competition 
is intensified when the number of element is increased. 
As a result, growth rate of real systems usually decreases 
with the increase of the system size, exhibiting the so- 
called density- dependent growth [26|. In fact, competi- 
tive dynamics has been found to dominate the evolution 
of variety of social 2?J , biological [28[ and economic [29j 
network systems. Thus, with respect to the AD growth 
rule, to be more realistic, P a should be a decreasing func- 
tion of the number of nodes in the network, rather than 
a constant [13] • In this paper, we introduce the Logis- 
tic population equation [26] into the AD growth rule. 
As a result, the probabilities of addition (deletion) be- 
comes a decreasing (increasing) function of the network 
size. Based on this dynamical AD growth rule, grow- 
ing networks models with and without PA are investi- 
gated, respectively. In the present models, networks ex- 
hibit the density-dependent growth and all the observed 
FTDDs are created. Moreover, it is found that the net- 
works exhibit nonstationary degree distributions, chang- 
ing from power-law to exponential one or from exponen- 
tial to Gaussian one. These results can be expected to 
shed some light on the formation and evolution of real 
complex real- world networks. 



II. NETWORK MODELS WITH LOGISTIC AD 
GROWTH RULE 

The Logistic equation (2(| was proposed by Pierre Ver- 
hulst for the analysis of population competitive dynam- 
ics. Verhulst assumed that during the population growth, 
as a result of the competition for the limited resources, 
the death rate per individual is a lineally increasing func- 
tion of the number of population. Then the Logistic 
equation can be written as 



dN(t) 
dt 



[a-cN(t)]N(t) =a[l 



N(t) 
K 



Mt) , (1) 



where N(t) denotes the number of population at time i; 
a and c are constants; a is the birth rate per individual 
and cN(t) stands for the death rate per individual at 
time t; K = a/c is called carrying capacity ,3d, [3l| which 
represents the largest population the environment can 
support due to limited resources. 

A more realistic description to the real-network's evo- 
lution can be achieved by the introduction of such Lo- 
gistic dynamics into the AD growth rule. To do this, we 
obtain the probability of addition P a and the probability 
of deletion Pd by the following two equations 



Pa 
P, 



aN(t) 
cN(t) 2 



(2) 



and 



which yield 



Pa+Pd = 1, 



and 



Pa(t) = 



Pd{t) 



K 



+ cN{t) K + N{t) 



N(t) 



cN{t) 

+ cN(t) ~ K + N{t) ' 



(3) 



(4) 



(5) 



Where N(t) stands for the number of nodes at time t and 
K is the parameter which denotes the maximum nodes 
in the network. 

Then the Logistic AD model can be defined as follows: 
We start from to isolated nodes, which act as the nu- 
clear of a growing network. At each time step, either a 
new node is added to the network with probability P a or 
a randomly selected old node is removed from the net- 
work with probability Pd = 1 — P a , where P a and Pd are 
determined by Eq. Q and Eq. (O, respectively. 

When a new node is added to the network, there are 
still two ways for it to attach to the existing nodes in the 
network. One way is to randomly choose m nodes to set 
up connections (growing network without PA) [TP . fl2j , 
and the other way is to preferentially select m nodes to 
connect by means of preferential probability in the BA 
model [H, which reads 



n i; = 



1 



i 



(6) 



where ki is the degree of the ith node (growing network 
with PA). Here, we should note that in order to give 
chance for isolated nodes to receive a new edge we choose 
Hi proportional to ki + 1 [H| . 

One can find from Eq. ^ and Eq. (|5|) that when 
N(t) < K, P a (t) > P d (t). This means that the networks 
grow rapidly at the initial stages of their evolution. In 
fact, in a Logistic model [jt| : 



dN(t) 
dt 



K - N(t) 



when N(t) < K, 



the solution is: 



dN(t) 
dt 



K 



aN(t), 



N(t), 



N (t) = C exp (at) , 



(7) 



(8) 



(9) 



where Co is a integral constant. Indeed, this kind of 
exponential growth has been observed in many newly 
emerged real- world networks, such as the World-Wide- 
Web (WWW) and the Internet [H,[20|. This indicates 
that the rapid growth of some real systems in their young 
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age is well described by the Logistic AD model. As N(t) 
increase, P a (t) decrease while Pd(t) increase, the growth 
of network is slowed down. In the limit of large t, 



and 



lim P a (t) = lim P d (t) = -. 



lim N (£) = K, 



(10) 



(11) 



i.e., the network reach a steady state and the number 
of nodes has the upper limit K. The above analysis 
imply that networks in the Logistic AD model exhibit 
density-dependent growth characterizing the evolution of 
many real systems, not captured in most previous net- 
work models. 

We investigated the degree distribution of the networks 
by extensive computer simulations. In the simulation, 
we set mo = m = 5 and K = 100000. Network is 
left to evolve until the steady state is reached. Cumu- 
lative degree distributions of growing networks with and 
without PA at different time step are given in Fig. [TJa) 
and Fig. QJb), respectively. We found that in the Logis- 
tic AD model, the networks exhibit nonstationary degree 
distribution before the steady state is reached. This is 
illustrated in Fig. [TJ For the growing network with PA, 
as time goes on, P(k) of the network undergos a pro- 
cess of transition, which can be roughly separated into 
several stages: (1) At the earlier stage of network evo- 
lution, i.e., when t < 72000, the network exhibits vari- 
ous PLDDs with different power-law exponents. In ad- 
dition, the exponent increases with time. Particularly, 
in the asymptotic cases of P a being close to 1, the net- 
work is almost equivalent to the well-known BA model, 
thus it have PLDD with power-law exponents 7 = 3. (2) 
72000 < t < 1000000. P(k) of the network is truncated 
by an exponential cutoff and the network exhibits the 
TPLDD during this stage. Generally, if a node needs 
to obtain the high degree, it must exists in network with 
enough long time. While, longer do the nodes live, higher 
is the probability that they are deleted. Moreover, with t 
increasing, Pd increases. Thus the surviving probability 
of the nodes with high degrees is greatly reduced, and the 
degree distribution evolves with time, gradually changing 
from a power-law form to a power-law with a exponen- 
tial cutoff, then to a exponential one. (3) In the limit of 
large t, e.g., when t > 1000000, the network shows a well 
shaped stationary EDD [see Fig. QJa)]. Speciously, this 
asymptotic case seems to be a non-growing network (or 
a very slowly growing one). The case of the preferential 
attachment on a non-growing network was considered in 
Ref. [l9[ where it was found that P(k) is not stationary, 
changing from a power-law type directly to is a Gaussian 
one. On the other hand, for the growing network without 
PA, as time goes on, P(k) of the network exhibits a con- 
tinuous transition from the EDD to a variety of TEDDs 
which prove to be a series of intermediate states, and in 
the limit of large t, e.g., when t > 1045700, the network 
exhibits a well shaped stationary GDD [see Fig. [ljb)]. 
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FIG. 1: (a) Cumulative degree distribution of the growing 
network with PA, for different time steps t, in logarithmic 
scales. The dash line is power-law fit and the solid line is 
exponential fit. (b) Cumulative degree distribution of the 
growing network without PA, for different time steps t, in 
semi-logarithmic scales. The dash line is exponential fit and 
the solid line is Gaussian fit. Results in both (a) and (b) are 
based on many independent realizations. 



The results of Fig. [T] reveal a nonstationary behavior 
for the degree distribution of real systems. During the 
evolution of these systems, various events take place on 
different timescales. These events include, for instance, 
the addition and deletion of nodes, the creating and 
rewiring of edges between existing nodes (internal edges), 
and so on. In general, the timescale on which a node 
join or leave the network (the addition and deletion of 
nodes) may be much longer than the timescale on which 
other events take place. For example, in a social net- 
work, creating internal edges means that the individuals 
make new friends, which happens on the timescale that 
can be as short as hours or days. While the timescale on 
which individuals are born or die is typically some years 
or decades. So that our model, which is based on the 
mechanism of addition and deletion, gives the prediction 
of a long-run behavior of real systems. For this reason, 
this kind of nonstationary FTDD of real systems, which 
can be difficult to observe in a relatively short time inter- 
val, has been neglected in most previous studies [l], [p| . 
However, in recent years, the Internet and the WWW is 
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in the initial and rapid growth stage of their evolutions. 
Perhaps, the observations for their degree distributions 
might provide some clues for our results. The WWW at 
the document level had grown at least five times larger 
during the two years time delay between the first and 
last web crawl, and the degree distributions obtained at 
different time are reported to be of a power-law form 
[32l - l34l |. The power- law exponent, though seems to be in- 
variable for the in-degree distributions, has an increasing 
tendency with the sample size or time for the out-degree 
distribution, changing from 2.45 in 1999 to 2.72 in 2000 
[32l - l34j ]. Another clue is refereed to the Internet on the 
inter-domain level. The degree distribution of the Inter- 
net is also of power-law form. The exponent, however, 
suffer a little change from 2.15 in November 1997, to 2.16 
in April 1998, then to 2.20 in December 1998 [H]. In fact, 
it is hard to achieve this precision. One may estimate the 
value of the exponent using the highest degrees and the 
system's size 36]. Such estimations also confirm the re- 
ported values. For November 1977, we gets 7 ~ 2.22, for 
April 1998, 7 ~ 2.24, and for December 1998, 7 ~ 2.26. 
Of course, exist additional factors (nonlinear PA [36[ and 
accelerating growth of edges [37], etc.) that may change 
these values. Finally, we are glad to point out that in a 
recent work [14J, such nonstationary degree distribution 
was also reported when the creation of internal edges was 
considered in growing and static networks. 

It is well known that the structure of complex net- 
works has strong effects on their function. Therefor such 
neglect, at least in some case, may be dangerous because 
networks with different degree distribution act differently 
if we consider the dynamic processes taking place on 
them. This is especially the case, for example, in the 
strategy designing for the prevention of virus spreading 
or for the defense of network attacking |38l - l40l |. Pro- 
vided that, after some years, the WWW has gradually 
changed into an exponential network, but is still treated 
as a scale-free one, what will happen to this network of 



high technological importance if it runs under the an- 
tivirus strategy designed for a scale-free network? 



III. CONCLUSION 

In summary, we have proposed a dynamic growth rule 
which comprises the addition and deletion of nodes in 
the network. Based on this rule, by adding a node 
with a probability P a or deleting a node with probabil- 
ity Pd — 1 — P a at each time step, topological properties 
of growing networks with and without PA are studied 
respectively. The probabilities of addition and deletion 
are assumed to be determined by the Logistic population 
equation. In our model, networks exhibit the density- 
dependent growth which characterizes the evolution of 
various real systems, and all the observed degree distri- 
butions of real systems are created. Moreover, the net- 
works exhibit nonstationary degree distributions, chang- 
ing from power-law to exponential one or from exponen- 
tial to Gaussian one. It is discussed that this kind of non- 
stationary behavior can be unconspicuous in real-world 
networks. The result indicates that the degree distri- 
bution of real systems, which has been believed to be 
stationary in most previous studies, may in fact be non- 
stationary. 
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